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C/^ ; ABSTRACT 

■ Context. To investigate the earliest phases of star formation and study how newly-born stars interact with the surrounding medium, 

' we performed a line and continuum survey at NIR and mm-wavelengths of a sample of relatively isolated Bok globules. 

I Aims. We present a follow-up observational program of a star-forming site in the globule CB230. From narrow-band continuum 

observations of this site, we had discovered a bright [Fell] jet, which originates in the low-mass YSO CB230-A. We aim to investigate 
I ■ the physical properties of the region from where the jet is launched. 

O ' Methods. Our analysis was carried out using low-resolution NIR spectra acquired with the camera NICS at the TNG telescope, with 

f~l I JH and HK grisms and a 1 arcsec-wide slit. These observational data were complemented with infrared photometric data from the 

C/3 , Spitzer space telescope archive. 

Results. The relevant physical properties of CB230-A were constrained by SED fitting of fluxes from the NIR to the mm. The YSO 
' spectrum exhibits a significant number of atomic and molecular emission lines and absorption features. The characteristics of this 

^ ' spectrum suggest that we are observing a region in the close vicinity of CB230-A, i. e. its photosphere and/or an active accretion 

, disk. The spectra of the knots in the jet contain a large number of emission lines, including a rich set of [Fell] lines. Emission due to 

, . H2 and [Fell] are not spatially correlated, confirming that [Fell] and H2 are excited by dilferent mechanisms, in agreement with the 

models where [Fell] traces dissociative J-shocks and molecular hydrogen traces slower C-shocks. By using intensity ratios involving 
density-sensitive [Fell] lines, we estimated the electron densities along the jet to be 6 x 10^^-1 x 10"* cm"'. This indicates either high 
density post-shock regions of ionised gas or regions with a high degree of ionisation. 
I ■ Conclusions. By combining the present data with previously obtained maps at NIR- and mm-wavelengths, the emerging scenario is 

0^ ' that CB230-A is a Class O/I YSO driving an atomic jet that is observed to be almost monopolar probably due to its inclination to the 

' plane of the sky and the resulting higher extinction of its red side. This primary jet appears to be sufficiently energetic to open the 

00 , cavity visible in the NIR images and drive the large-scale molecular outflow observed at mm-wavelengths. CB230-A was revealed 

. to be a good location to test the innermost structure of accreting low-mass protostars. 

^ . Key words. Stars: formation - Stars: winds, outflows - ISM: jets and outflows - ISM: molecules - ISM: individual objects: CB230 

1- Introduction good example, given the simultaneous presence of almost all 

the ingredients expected in a typical star formation scenario: 

High-velocity jets driven by collimated winds generated by ysOs, a hot atomic jet, and a colder molecular outflow. One of 

Young Stellar Objects (YSOs) interact strongly with the sur- ^j^^ ySOs in the cloud (hereafter called CB230-A) was detected 

rounding medium, cleaning up the high-density material hosting jj^^ Near-Infrared (NIR) as well as at sub-millimetre wave- 

the star-forming process. Jets are also thought to remove angu- j^^^ ^yun & Clemens ESU fT994il Launhardt & Henning 

lar momentum from the accreting matter allowing lower angular 1997. Huard et al. 19991 Young et al. EOOl. It is also associ- 

momentum gas to carry on the build-up process of a new star ^^^^ ^^^^ ^ ^^ip^j^ ^^jA^^^ observed in CO (Yun & Clemens 

In addition, when they travel into the ambient cloud, jets create 1994^^ Using N2H+ observations, Chen et al. dSOOTj found a 

shocks that heat (up to thousands of K) and compress the gas, ^.j^^^ ^^^^^-^^ ^^^^-^^^ ^^^^^^ molecular core of ~ 8.8 km s ' 

and drive bipolar molecular outflows on larger spatial scales con- p^-i ^ ^^^^ direction as the line connecting CB230-A with 

taming colder (10-100 K) swept-up material. As a consequence, ^ ^^^^^^ ^ ^q„ ^^^^ j^^-^ j^^^^^ ^^^^^^ ^j^^^ CB230-A was 

the study of the acceleration and colhmation of jets is fundamen- ^^^^^^^^ 7 ^-^^^ ISOCAM and coincides with a possibly 

tal to understanding star formation. ^^^^^j^ ^JR source (CB230-B,C; see Fig.[B. 

Some of the most accessible places to survey jets are rep- 
resented by Bok globules (Bok & Reilly 1947), which are rel- The distance to CB230 was discussed by Launhardt & 
atively isolated molecular clouds mainly associated with low- Henning ( 119971 ). who placed it at 450 pc. This was derived 
mass star formation (e.g. Huard et al. 1999 and references by associating CB230 with a larger molecular cloud complex 
therein). The globule CB230 (Clemens & Barvainis |1988l ) is a (Cepheus Flare), based on their respective radial velocities. The 

authors estimated that the given distance is accurate to within 

Send offprint requests to: F. Massi (fmassi@arcetri.astro.it) 30%. Kun 01998( 1 discussed the distance to the Cepheus Flare 
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molecular cloud complex and for L1177 (the dark cloud hosting 
CB230) quotes a distance in the range 300-560 pc, where the 
lower value is the kinematic distance and the higher one is the 
distance to the nearby (but probably unrelated) reflection neb- 
ula vdB141. Although kinematic distances lower than 1 kpc are 
of doubtful reliability, the lowest value agrees with the distance 
(300 + 30 pc) derived for the dark clouds in the area by studying 
the cumulative distribution of stellar distance moduli (Wolf dia- 
gram). Young et al. ( I2006I I assumed d = 288 pc taking this value 
from Straizys et al. ( |1992| l. who derived the distance to a group 
of dark clouds (close to LI 177, but not including it) using pho- 
tometry of 79 stars. In particular, LI 167/Ll 174 were found to lie 
atd - 288 + 25 pc, in agreement with the lower value indicated 
by Kun d 19981 1 for LI 177. In summary, it seems that d - 450 pc 
(+ 30%) is a robust determination and 288 pc may be assumed 
as a safe lower limit. 

CB230-A is classified as a Class 0/Class I source (e. g., 
Froebrich i2005 ) . Estimates of the bolometric luminosity, based 
on the IRAS fluxes, range from 7.7 Lq (Froebrich |2005l l to 
12 L0 (Launhardt & Henning '1997), both for a distance of 
450 pc, i. e. the distance we adopt. From our sub-mm obser- 
vations with SCUBA (Brand et al. 1200^ . we derive a gas mass 
of 3.4 M0, which agrees within a factor of 2 with the mass es- 
timated by Launhardt ( 20011 1 from observations at 1.3 mm with 
flie IRAM 30-m. Launhardt (2001) also found a disk of mass 
~ 0.1 M0 around the protostar, based on observations at 1mm 
continuum with the Owens Valley Radio Interferometer. By us- 
ing three different evolutionary models, Froebrich (|2005l l de- 
rived an age in the range 1.7 x 10"^-2.2 x 10^ yrs and a final 
stellar mass in the range 0.3-0.9 Mq. 

In the framework of our ongoing multi-frequency study of 
star formation in Bok globules (Massi et al. 2004; Codella et 
al. I2006i l. we carried out narrow-band NIR observations cen- 
tred on the 1.644 fim [Fell] and 2.122 jjm H2 lines. As a re- 
sult, we detected a bright [Fell] jet in CB230, which (i) orig- 
inates in CB230-A and extends for ~ 0.02 pc, (ii) is directed 
along the same direction as the molecular outflows axis (N-S), 
and (iii) is defined by two knots (called kl and k2) superim- 
posed on fainter elongated emission observed also in H2 (Massi 
et al. [20041 1. Figure [1] illustrates this by displaying an overlay 
of the [Fell] image and the CO(3-2) bipolar outflow as mapped 
by Brand et al. ( |2008] l. These [Fell] and H2 NIR fines are par- 
ticularly useful, since [Fell] emission traces fast and dissociative 
J-shocks, which should outline the inner jet-channel, whereas H2 
emission is expected to trace slower C-shocks that probably orig- 
inate in bow-shocks, probing the region of interaction between 
the jet and the ambient gas. The CB230 NIR maps therefore sup- 
port a complex scenario, which, to be confirmed, requires further 
[Fell] and H2 spectroscopic NIR observations; these data were 
obtained and enabled us to constrain the physical properties of 
the gas (such as extinction and excitation temperature) in the re- 
gion where the jet interacts with the surrounding dense gas close 
to its origin, i. e. close to the central system, which consists of a 
protostar and, presumably, an accretion disk. 

2. Observations and data reduction 

The spectroscopic observations were carried out on the night of 
August 3, 2004, with the NIR camera NICS at the Telescopio 
Nazionale Galileo (TNG) at La Palma. We obtained spectra 
through the grisms 7// ( 1 . 1 5 - 1 .75 //m) and //T*: ( 1 .40-2 .50 fim), 
using a l"-wide slit. This provided a spectral resolution of 
R~500. The slit was aligned along the [Fell] jet defined by the 
kl and k2 knots, hence encompassing also the driving source 




right ascension (JSOOO) 

Fig. 1. Continuum-subtracted [Fell] map (grey scale; Massi et 
al. 2004) of CB230: the white contours mark the brightest emis- 
sion indicating the YSO CB230-A. The label CB230-B,C mai'k 
two NIR objects, that could make a triple system with CB230-A. 
The labels kl and k2 indicate the two [Fell] knots defining the 
jet. Also visible is an extended diffuse component showing the 
cavity opened by the YSO mass ejection. In contours, the CO(3- 
2) integrated emission of the blue (continuous) and red (dashed) 
lobes of the molecular outflow driven by CB230-A (Brand et 
al. 2008; JCMT-data). The integration intervals are from -6 to 
-1-1.37 km s"' (blue lobe) and from -1-4.68 to H-14 km s"^ (red 
lobe). Contours range from 3 to 20 K km s ' in steps of 1 K km 
s '. The blue lobe extends for ~ 120", the red one is roughly 
half as long. 



(the YSO). Six spectra of the target were obtained through the 
grism JH in an ABBAAB set of cycles (i. e. integrating when 
the source is in a position A on the detector, then nodding the 
telescope along the slit projection on the sky to move the source 
to a different position B on the detector, integrating again and 
reversing the cycle, etc.). The integration time for each spectrum 
was 900 s. Two spectra were taken through the grism HK in a 
single AB cycle, with integration times of 600 s each. The star 
HIP97033 (spectral type GO V) was observed in between the 
JH- and the HK-cycles on the source, by acquiring four spectra 
in an ABBA cycle with 2 s of integration each, through both 
the JH- and the HK-gnsm. Other spectra of this star and of 
HIP108772 (spectral type 09 V) were taken at the beginning of 
the night, but they could not be used to correct the target spec- 
tra for telluric absorption because they were observed at least 5 
hours before the target. However, they were used to obtain the 
intrinsic spectrum of HIP97033 and perform other tests on the 
telluric correction. The weather conditions did not remain stable 
and observations of the target could be carried out only in the 
late night due to clouds. 
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Data reduction was done using IRAlQ All frames were flat- 
field corrected by using spectra of a halogen lamp, following 
the prescriptions given on the TNG web page, i. e. the halogen 
frames were normalised row -by-row by dividing by a low-degree 
polynomial fit of the background. All frames were wavelength- 
calibrated and rectified by using reference spectra of an Argon 
lamp. Then, each pair of subsequent frames (i. e. the A and B ex- 
posures in AB- or BA-cycles) were subtracted to remove back- 
ground emission and other biases. Since a single frame (i. e. A 
or B in AB- or BA-cycles) is an image where the source spec- 
trum appears as a row of bright pixels, after subtraction a frame 
is obtained that contains two spectra of the source, one posi- 
tive and one negative with respect to the average background 
level. By careful examination of the [Fell] images in Massi et al. 
( 1200 4'). we estimated the extent of the two emission knots and 
the YSO in the spatial direction, and extracted their spectra from 
each subtracted frame (two per frame). Two spectra of HIP97033 
per subtracted frame were also extracted. 

Using the solar spectrum as an input, a synthetic spectrum 
was constructed (following Maiolino et al. 119961 ) that repre- 
sented an arbitrarily-scaled version of the intrinsic spectrum 
of HIP97033 (assuming a radial velocity of -21.6 km s"' and 
a rotational velocity of 10 km s ', as listed in the SIMBAD 
database). For all extracted spectra of HIP97033, the wavelength 
calibration was then further refined by aligning the intrinsic ab- 
sorption lines in each stellar spectrum with the corresponding 
ones in the synthetic (reference) spectrum. This is required be- 
cause moving the telescope produces small shifts in the wave- 
length scale of each frame (i. e. there exists a shift between 
the wavelength scale in the Ar-lamp frames and that in each 
of the source frames). Finally, all other spectra were further 
wavelength-corrected (to minimize this instrumental effect) by 
shifting their wavelength scale to align their atmospheric absorp- 
tion lines with those in the final stellar spectra. This correction 
translated into shifts of ~ 10"^ jum (i. e. ~ 10 A, independent 
of wavelength) both for the telluric standard and the targets. All 
spectra of a given source were then coadded to obtain the high- 
est signal-to-noise ratio. The atmospheric transmission was cor- 
rected for as explained in Maiolino et al. ( 1996 ), i. e., by dividing 
the spectra of the knots and the YSO by that of the telluric stan- 
dard and then multiplying them by the synthetic intrinsic spec- 
trum of the star obtained from the solar spectrum. 

2.1. Flux calibration 

We first dereddened the spectra over the entire wavelength range 
by adopting the reddening law of Rieke & Lebofsky ( 1985]l. We 
tentatively derived Ay from [Fell] line ratios (see Sect. 13.11 the 
uncertainties affecting the obtained values are further discussed 
in Sect. 14.31 1. Once Ay was known, the reddening correction was 
obtained by dividing each spectrum by a curve proportional to 
the extinction as a function of the wavelength. Then, the dered- 
dened spectra were flux-calibrated according to the following 
prescriptions. For the knots (i. e. Figs. |2] and [3]), we required that 
the total flux of the [Fell] line at 1.645 jjm equalled the narrow- 
band photometric determinations of Massi et al. ([20041) . 

For CB230-A (i. e. Fig.H]), each of the J, H, and K segments 
of the dereddened YSO spectrum was first multiplied by the to- 
tal transmission of the 2MASS system (atmosphere excluded); 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



the derived flux densities were then integrated over each of the 
J, H, and K passbands and divided by the integral of the total 
2MASS transmission over the corresponding passband. The re- 
sulting mean flux densities were eventually set to be equal to the 
flux densities derived from the 2MASS magnitudes of the YSO 
(2MASS21173862+6817340,y = 14.085 + 0.062,//= 11.712, 
K = 10.487) after being dereddened using the same Ay used to 
deredden the entire spectrum. We then obtained two conversion 
factors (counts to flux density) from JH (one from J and one 
from H) and two conversion factors from HK (one from H and 
one from K). 

Obviously, the two conversion factors for each of the two 
YSO spectra (i. e. JH and HK) should coincide. For the HK 
spectrum, the conversion factors obtained from H and K agree 
to within a few percent, whereas for the JH spectrum, those ob- 
tained from J and H differ by ~ 30%. The two different calibra- 
tion factors within JH probably reflect errors in the photometry 
due to the presence of diffuse emission towards the YSO and 
possible time-variability slightly affecting the SED (spectra and 
2MASS data were obviously not taken simultaneously). 

However, the calibrated YSO spectra have different contin- 
uum emission slopes within the overlapping H interval, with dif- 
ferences in flux of up to ~ 20 %. The same problem is found 
when calibrating with the integrated flux in the [Fell] line at 
1.645 /urn, which produces differences in the slope of the under- 
lying continuum towards the knots of up to ~ 30 %. Inaccuracies 
in the telluric correction are discussed in AppendixlAland do not 
appear to affect the spectra significantly. The most likely cause 
is the intense, extended emission around the source, which could 
have prevented an accurate placement of the slit on the central 
source. 

2.2. Spitzer observations 

Spitzer observations of CB230 were taken from the Spitzer pub- 
lic archive. They consisted of IRAC images at the four wave- 
lengths (3.6, 4.5, 5.8, and 8 fim) and MIPS images at 24 and 
70 jum. The observations were part of the GTO program 124, 
"IRAC and MIPS Imaging and IRS Spectroscopy of Pre and 
Post Main Sequence Stellar Systems" (Principal Investigator G. 
Gehrz). Overviews of IRAC and MIPS were given by Fazio et 
al. (12004 1) and Rieke et al. (12004), respectively. 

We retrieved all post-BCD (Basic Calibrated Data) images, 
which were produced by the pipeline version S 14.0.0 (IRAC), 
S 16. 1.0 (MIPS at 24 yum), and S 16.0.1 (MIPS at 70 yum). All 
MIPS data were obtained in the photometric small field modes 
(fine scale at 70 //m). After checking all images for possible 
problems as described in the instrument Data Handbooks 0, we 
decided not to repeat the reduction steps starting from the BCD 
images, and photometry was carried out on the post-BCD im- 
ages. 

CB230-A and CB230-B(,C) are resolved in the IRAC and 
MIPS 24 jum images, but blended in the MIPS 70 fim im- 
ages. CB230-B and CB230-C are unresolved even in the IRAC 
images, so they were considered to be a single source (here- 
after, CB230-B-HC). Aperture photometry of both CB230-A and 
CB230-B-1-C was performed using DAOPHOT in IRAF We se- 
lected aperture radii of 2 pixels (~ 2" [IRAC], 3" [MIPS 24 //m], 
and 35" [MIPS 70 //m]). The background contribution was mea- 
sured within annuli with inner radius of 2 pixels and width of 
4 pixels (IRAC), inner radius of 6" and width of 7" (MIPS at 
24 fim), and inner radius of 39" and width of 26" (MIPS at 

" http://ssc.spitzer.caltech.edu/ 
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Table 1. Spitzer magnitudes of CB230-A and CB230-B+C. 



Band 


CB230-A 


CB230-B+C 


(jum) 


(mag) 


(mag) 


3.6 


9.101 +0.004 


12.62 ± 0.04 


4.5 


8.537 + 0.006 


12.06 ± 0.04 


5.8 


8.033 + 0.006 


11.49 ±0.06 


8 


7.88 ±0.01 


10.21 ±0.01 


24 


2.847 ± 0.004 


4.48 ± 0.02 


70° 


-3.75 


±0.07 



" CB230-A and CB230-B+C are spatially unresolved in this band and 
appear as a single source. 

24 jum). We adopted aperture corrections given on the Spitzer 
web page0, for the used apertures and sky annuli. In the case of 
the MIPS 24 fim image, the aperture photometry for each of the 
two sources was performed after subtraction of the other source 
by PSF fitting. The photometric uncertainties were computed as 
described by Reach et al. ( |20051 l and Gordon et al. jlOOlj . We 
note that due to flux density nonlinearities, aperture photometry 
is expected to underestimate the MIPS 70 jum flux, probably by 
~ 20 % (Gordon et al. |2007| l. The fluxes were converted into 
magnitudes in the IRAC (Reach et aI. I2005l l and MIPS (see the 
MIPS handbook) systems and are listed in Table[T] 

3. Results 

3. 1 . Emission lines 

The complete JHK spectra of the knots kl and k2 and of the 
YSO CB230-A are shown in Figs.|2l|3] andg] respectively. The 
spectra were dereddened using the visual extinctions derived 
from [Fell] line ratios (see below). In Table |2] we list fluxes 
and identification for all the detected emission lines at the three 
positions (kl, k2, and CB230-A). We note that we considered 
detections to be only those lines with an S/N > 3 and with line 
widths of comparable or larger than the instrumental resolution. 
Table |2] also contains a number of tentative detections (also la- 
belled in Figs.|2]|3] and|4]l with a S/N ^ 2. The three spectra show 
a rich set of [Fell] emission lines, as well as the H2 emission 
line at 2.12 fj.m: in particular, the knot kl exhibits a large num- 
ber of bright emission features. These confirm the presence of 
H2 emission and unveil the structure of the [Fell] line-emission, 
providing improvements to the interpretations of narrow-band 
filter data by Massi et al. ( 120041 ). 

3.2. YSO spectral sliape 

All spectra show an underlying continuum emission, whose 
shape (before correction for extinction) is similar in all direc- 
tions, i. e. towards both the YSO and the knots. However, in 
the YSO spectrum, the continuum emission is more intense and 
exhibits several emission and absorption features. Longward of 
1.3 fim, the beginning of an absorption band can be seen in 
all spectra, which can be identified as being due to water (see 
Sect. |4.4| l. A check of the spectra uncorrected for telluric absorp- 
tion confirms that this band is intrinsic and not a telluric feature. 
The continuum emission from the knots most likely originates 
from the scattering of radiation by dust grains. The ubiquitous 
presence of this water absorption band indicates that the illumi- 
nation source is the YSO. Another absorption band clearly vis- 
ible longward of 2.29 yum, is easily identifiable as the CO(2,0) 

' http://ssc.spitzer.caltech.edu/ 



band. This clearly indicates that the central star photosphere 
and/or the innermost region of the circumstellar disk are being 
observed through a cavity produced by the outflow. 

4. The YSO 

4.1. YSO classification 

The classification of both CB230-A and CB230-B+C can be 
checked by using Spitzer data to derive their spectral indices be- 
tween 2 and 24 fim. Assuming that AF^ oc A", we estimated that 
a ~2 (CB230-A) and a ~ 2.4 (CB230-B-HC) from the 2MASS 
K fluxes and MIPS 24 jum fluxes. Both objects can be then iden- 
tified to be Class I sources. 

The Mid-Infrared (MIR) colours also, derived from the mag- 
nitudes given in Table [1] can provide information about the evo- 
lutionary stage of the two objects. In a [3.6]-[4.5] versus [5.8]- 
[8.0] colour-colour plot (CCP), both lie outside the "disk do- 
main" (i. e. the plot region dominated by Class II sources) from 
AUen et al. ( |2004| l. CB230-A colours would be consistent both 
with those of a reddened photosphere (although for an Ay of 
at least 40 mag), as well as with those of an accreting proto- 
star ("stage" I source, according to the classification proposed 
by Robitaille et al. |2006] l. CB230-BH-C colours are located in a 
region occupied mainly by both stage II sources (in practise, they 
are mostly Class II sources) and stage I sources. The information 
from a [3.6]-[5.8] versus [8.0]-[24.0] CCP is more helpful. Both 
CB230-A and CB230-B-I-C colours are consistent with those of 
stage I sources or heavily reddened {Ay > 20 - 30 mag) stage II 
sources. 

4.2. YSOSED 

As a further step towards an understanding of the evolutionary 
stage and physical properties of CB230-A, we used a more com- 
plete set of observational data spanning the wavelength interval 
from the NIR to the mm. Model SEDs could then be fitted to 
the data and the relevant physical properties of the YSO derived. 
The source fluxes at different wavelengths, and the effective or 
estimated aperture radii for which they were measured, are listed 
in Table |3] 

By integrating the SED, we found that the ratio of the lu- 
minosity for wavelengths longer than 350 fim to the bolometric 
luminosity is ~ 0.01, which would lead to its classification as a 
Class source (Andre et al. 2000). Given also its properties in 
the NIR to MIR, we can therefore understand its classification as 
a transition Class 0/Class I source. 

To fit the SED o f CB2 30-A, we used the onHne fitting tool 
by Robitaille et al. ( |2007l ). This exploits a grid of 20000 pro- 
tostar models for which the radiation transfer has been solved, 
accounting for the central star, envelope and disk and includ- 
ing cavities and scattering. A wide range of physical conditions 
in each component is covered (Robitaille et al. 2007). We con- 
strained the distance to be in the range 280-600 pc (see Sect.lTJ. 
All models (i. e., 16) within a range ~ xlest ^ ^P^^ 
apoint) are shown in Fig.|5j where xl^^^ - 233.47 is the value 
(per datapoint) from the best fit. Clearly, these appear to fit the 
SED reasonably on a by-eye basis. The high values ofx^ are due 
partly to having underestimated the uncertainties (in most bands, 
we adopted the photometric errors only). In fact, we found that 
by increasing the errors in the IRAC and 1.3 mm fluxes to 10 
%, the of the best-fit model decreased to 315 (i. e. 20 per 
datapoint), whereas the range of physical parameters is not sig- 
nificantly changed. Nevertheless, we also checked that none of 
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Table 2. Emission lines observed towards the knots kl and k2 and the YSO CB230-A: detections and tentative detections. The 
values between brackets are the measured rm.s. 



kl" k2" CB230-A' 



Transition /l* (jum) F (10 erg cm ^ s ') 



rinoTTl o-^ Pi 

[rSilJ a -L'7/2 ^ 


a 




1 T^A7 
i .ZJU / 


no (W 
/b.O(,Z.U) 


lzj.(5(4.(5; 




[reilj a 1^1/2 — 


a 


JJl/2 


i.Z/U4 






/.y(i.D; 


ri"' TTT 4 1 — V 

[Fell] a^ D3/2 - 


6 

a 


D3/2 


1.2788 


10.8(2.0) 




14.2(1.6) 


[Fell] D5/2 - 


a* 


D5/2 


1.2943 


11.1(2.9) 




7.7(2.0) 


[Fell] D7/2 - 


a* 


JJ7/2 


1 QonA 
l.j/UD 




1 A/^ 1 \ 


lO.U(l.D) 


[Fell] a" D5/7 - 


a« 


D3/2 


1.3278 


7.6(2.3) 






[FeII]a2l„/2- 




H9/2 


1.5246 


5.0(1.3) 


10.7(5.2)'' 




[Fell] a" D5/2 - 


a* 


F9/2 


1.5335 


12.3(2.0) 


11.2(5.2)'' 


21.5(5.9) 


[Fell] a^ D3/2 - 


a'* 


F7/2 


1.5995 


5.7(1.7) 






[Fell] a-* D7/2 - 


a'* 


F9/2 


1.6436 


58.4(1.1) 


60.8(2.8) 


46.5(6.1) 


[Fell] a" Di/2 - 


a-* 


F5/2 


1.6638 


5.2(1.7) 






[Fell] a^ D5/2 - 


a"* 


F7/2 


1.6769'- 


4.7(1.1) 




10.9(6.3)" 


HI 11-4 






I.68I1'' 






15.7(6.3)" 


H2 1-0 S(l) 




2.122 


11.7(3.1) 


15.4(2.5) 


38.3(2.9) 


HI 7-4 






2.1661 






10.8(2.9) 



" The JHK spectra have been dereddened using the visual extinctions derived from the [Fell] 1.645/1. 321 line ratio (see Tab.[6](. '' The wavelengths 
of the [Fell] lines are taken from the NIST atomic spectra database. Blended with a telluric line. S/N ^ 2: tentative detection. " The HI and 
[Fell] lines at 1.68 jim are blended. 




1 ,25 13 1^59 t.62 1.65 1,68 1,71 

Fig. 2. (a) NIR {JHK) dereddened spectra (A,, = 14 mag) obtained with the TNG towards the [Fell] knot kl (see Fig.[T]i. Telluric 
lines have been marked, too. (b) and (c) Zoom-in on the spectral regions around the bright [Fell] emission lines at 1 .26 and 1 .64 jum, 
where several other weaker [Fell] lines have been detected. Table|2]provides the complete list of detections and tentative detections 
of emission lines. 
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1.2& 



1.3 
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Fig. 3. (a) NIR (JHK) dereddened spectra (Ay = 18.5 mag) obtained with the TNG telescope towards the [Fell] knot k2 (see Fig.[T]l. 
Telluric lines have been marked, too. (b) and (c) Zoom-in on the spectral regions where the bright [Fell] emission line at 1.26 jum 
and the H2 emission 2.12 fim line have been detected. Table |2] provides the complete list of detections and tentative detections of 
emission lines. 



the grid models is able to fit equally well the NIR-MIR and FIR- 
submm parts of the SED. In Table ID we list the most relevant 
physical parameters for the best-fit model and their range for the 
remaining 15 models within the above interval. 

Clearly, a picture emerges of a very young low-mass (~ 
0.5 Mq) central protostar with a high accretion rate. This agrees 
with its previous classification as a transition Class 0/Class I 
source. The age and mass determination also agree with the com- 
pilation by Froebrich (2005). Its inclination with respect to the 
line-of-sight is not large, as expected since only one side of a 
cavity is visible in the NIR. A circumstellar disk is present with 
a mass ~ 0.01 Mq. This is in accord with a disk mass ~ 0.1 M0 
inferred by Launhardt (2001) based on interferometric mm ob- 
servations, since Robitaille et al. (I2007I I checked that by using 
their fitting tool, the agreement between fitted disk masses and 
other disk mass determinations (when these are available), is 
always good to within one order of magnitude. However, we 
note that some of the best-fit models exhibit disk masses much 
smaller than ~ 0.01 M0, although the accretion rate is always 
quite high. This implies that even the available MIR data are un- 



able to constrain strongly the circumstellar disk physical proper- 
ties. 

Another problem is that all models appear to be underlumi- 
nous with respect to the estimated bolometric luminosity. This 
partly arises because of the range of fitting distances; after scal- 
ing all total luminosities instead to a distance of 450 pc (at which 
the reported bolometric luminosity estimates are derived) we 
achieve closer agreement between measurement and predictions 
of bolometric luminosity (see Table H]). Nevertheless, the total 
luminosity remains underestimated. This is clear from Fig. |5] 
where all model SEDs are always below the datapoints in the 
FIR-to-mm range. This reflects the intrinsic inadequacies of the 
model grid in describing the physical properties of CB230-A al- 
ready noted. These models do not account for example for the 
luminosity generated by the envelope accretion (Robitaille et al. 
2006), which may in fact be important in this case. 

We also attempted to fit the SED of CB230-B-HC with the 
same constraints as for CB230-A (but obviously without data- 
points in the FIR to mm range). The best-fit model has xl^^^ = 
17.84 (per datapoint) and its relevant physical parameters are 
listed in Table |5] along with those of the models (36) with 
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M 1.4 i,S i,S 2 M: 




1,2 1,25 1.3 1,6 1.64 1.66 

Fig. 4. (a) NIR (JHK) dereddened spectra (Ay - 8.4 mag) acquired by the TNG towards the YSO CB230-A, which drives the 
jet/outflow system (see Fig. [TJ. Telluric lines are indicated. Tables |2] and [T] provide the complete list of detections and tentative 
detections of emission and absorption lines, respectively. We note the absorption band at wavelengths less than 1 .2 fim, possibly due 
to VO, the three H2O absorption bands around 1.4, 1.9 and 2.4 yum, as well as the two CO absorption bands at 2.29 and 2.32 //m. 
Panels (b) and (c) present a zoom-in view of the spectral regions where the bright [Fell] emission lines at 1.26 and 1.64 fj.m are 
detected. Note that for the sake of clarity all the absorption features (detections and tentative detections), except for Mgl at 1 .50 yum 
and OH at 1.56 /im, stre labelled only in panels (b) and (c). 

- Xhest ^ ^P^^ datapoint). Unfortunately, the lack of mea- 
surements in the FIR and submm ranges prevents us from con- 
straining strongly many of the physical parameters. The best-fit 
model is that of a relatively evolved M~ 1 .5 Mq protostar with 
an edge-on disk that is no longer accreting. This can probably 
be discarded, based on the interferometric mm measurements by 
Launhardt (2001) that set an upper limit of 0.006 M© to the disk 
mass of CB230-B+C. However, there still remains a combina- 
tion of Class I and Class II sources with central masses in the 
range ~ 0.1 - 1.5 M© that are able to fit the SED of CB230- 
B-hC. Given the possible binary nature of this source, more data 
are also required to aiTive at a more definite identification. This 
source will not be discussed further in the following. 

4.3. Scattering and extinction in the NIR 

In principle, the NIR spectrum of the YSO provides a powerful 
means to constrain further its total photospheric absorption. The 



most intense emission lines are due to three [Fell] transitions at 
1.257 yum, 1.321 yum, and 1.645 yum. Since these [Fell] lines cor- 
respond to the same upper energy level (see Table |2]i, the line 
ratios 1.645/1.321 yum and 1.645/1.257 yum can be used to es- 
timate the extinction towards the [Fell] -emitting regions. The 
three lines are within the observed wavelength range of grism 
JH; determination of these ratios is therefore straightforward 
and no intercalibration that exploits partly overlapping spectral 
intervals is required. The fluxes were computed from Gaussian 
fits to the lines and the noise was measured in small adjoining 
spectral intervals. Finally, the extinction was derived towards 
the YSO and the two knots following Nisini et al. (2005b), and 
is listed in Table |6] The quoted uncertainties were obtained by 
propagating the noise. Clearly, there appear to be some discrep- 
ancies in the extinction estimated from the two line ratios. Nisini 
et al. (|2005bl l pointed out that these discrepancies are generally 
found when deriving the extinction from the two line ratios, such 
that the 1 .645/1 .257-ratio appears to overestimate Ay. They rec- 
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Table 3. Fluxes of CB230-A used for SED fitting. Those for 
CB230-B-I-C are also listed when the two (unresolved) sources 
are not blended with CB230-A. 



Band 


CB230-A 


CB230-B+C 


Aperture radius 


(fim) 


(mJy) 


(mJy) 


(") 


1.235" 


3.702 ±0.211 


< 0.748 


4 


1.662" 


21.159 


1.128 ±0.124 


4 


2.159" 


42.57 


3.893 ±0.218 


4 


3.6* 


64.3 ± 0.2 


2.51 ±0.08 


2.4 


4.5* 


69.1 ±0.4 


2.70 ± 0.09 


2.4 


5.8* 


70.4 ± 0.4 


2.9 ±0.2 


2.4 


8* 


45.4 ± 0.4 


5.27 ± 0.07 


2.4 


12" 


< 250 


67 


24'' 


521 ±2 


116±2 


3 


25" 


680: 


t70 


67 


60" 


11750 


±700 


95 


70"' 


24566 : 


t 1500 


35 


100" 


33530: 


t2000 


134 


450" 


10430 


50 


850" 


2330 


20 


1300^ 


221 


±5 


6 



" 2MASS catalogue; * Spitzer/IRAC, this wor k; " I RAS PSC; '' 
Spitzer/MIPS, this work; " SCUBA, Brand et al. ( GOOSl l; ^ IRAM 30- 
m, Launhardt & Henning I il997t ; 




1\ (.wnj 

Fig. 5. Model SEDs (full colour lines in the electronic edition) 
that best fit the fluxes of CB230-A (full dots; the full triangle 
marks an upper limit), within a range ~^best ^ ^P^^ data- 
point). 

The dashed line marks the stellar photosphere of the best-fit model. 



ommend using the 1. 645/1. 321-ratio, ehich they found provides 
results that are more consistent with other independent determi- 
nations. However, in our case it is evident that the two values 
of Ay are always equal within the uncertainties, that appear then 
dominated by the line signal-to-noise ratio, rather than by the 
uncertainty in the transition probabilities suggested by Nisini et 
al. ( |2005bl ). However, in Appendix IB] we propose an alternative 
explanation for the disagreement between the extinction values 
obtained from the two line ratios: the 1 .257 fim line intensity can 
be underestimated because of unaccounted-for telluric absorp- 
tion. In the following, we follow the prescription of Nisini et al. 
(I2005bl) and use extinctions derived from the 1.645/1. 321-ratio. 

The error in line ratios due to inaccurate correction of the at- 
mospheric absorption (and large-scale differences in system sen- 
sitivity) can be estimated to be ~ 10%, according to the discus- 
sion in AppendixlAl This would cause the extinction A i 544 to 



Table 4. Relevant physical parameters of the protostellar model 
that most accurately reproduces the SED of CB230-A and their 
ranges for the remaining 15 models in the interval "A'best ^ 
250 (per datapoint). 



Physical parameter 


Best fit 


Range 


x' 


3735.51 


3804.37-7668.63 


Interstellar extinction Ay (mag) 


0.00 


3.79-18.02 


Distance (pc) 


295 


282-513 


Age (yr) 


5.2 X 10" 


2.7 X 10*-W^ 


Central mass (Mq) 


0.4 


0.15-0.5 


Star radius (Rq) 


4 


3.4-4.9 


Star temperature (K) 


3600 


2908-3759 


Accretion rate (Mq yr"') 


10-* 


10-^-10-* 


Envelope outer radius (AU) 


5784 


1968-7757 


Cavity aperture angle (deg) 


23 


21-29 


Disk mass (Mq) 


0.03 


7.7 X 10-^-0.05 


Disk inner radius (AU) 


0.1 


0.06-0.3 


Disk accretion rate (10^* Mq yr"') 


2.3 


0.0077-3.6 


Inclination to the line of sight (deg) 


31.8 


18.2-31.8 


YSO inner extinction Ay" (mag) 


48.09 


0.09-37.33 


Total photospheric extinction 






Av+A'y (mag) 


48.09 


6.5-42 


Total system luminosity (Lq) 


2.585 


0.846-4.287 


Total system luminosity 






scaled to d = 450 pc (Lq) 


6 


2.2-5.7 


I. e., the extinction from the outer edj 


;e of the circumstellar envelope 



to the stellar surface along the line of sight. 

Table 5. Relevant physical parameters of the protostellar model 
that most accurately reproduces the SED of CB230-B+C and 
their ranges for the remaining 36 models in the interval - 
-Vbest < 30 (per datapoint). 



Physical parameter 


Best fit 


Range 


x' 


142.72 


242.36-382.14 


Interstellar extinction Ay (mag) 


15.04 


5.26-24.96 


Distance (pc) 


525 


316-602 


Age (yr) 


1.9 X 10'' 


5.2 X 10^-10'' 


Central mass (Mq) 


1.4 


0.10-1.59 


Star radius (R©) 


2.4 


1.6^.3 


Star temperature (K) 


4540 


2559-5655 


Accretion rate (Mq yr"') 





0-4.5 X 10-"' 


Disk mass (Mq) 


0.037 


2.7 X 10-^.029 


Disk inner radius (AU) 


36 


0.044-16.6 


Disk accretion rate (10"' Mq yr"') 


4.0 


0.00019-1.69 


Inclination to the line of sight (deg) 


87.13 


18.2-87.13 


YSO inner extinction Ay" (mag) 
Total system luminosity (Lq) 


459.4 


4x 10^^-182 


8.230 


0.33-6.91 



" I. e., the extinction from the outer edge of the circumstellar envelope 
to the stellar surface along the line of sight. 



be over- or under-estimated by ~ 0.2 mag, which translates into 
AAy ~ 1.2 mag. Comparing this with the values listed in Table|6] 
shows that in all cases but one (i. e. kl), the uncertainty in the 
reddening due to the signal-to-noise ratio of the single lines is 
larger or of the order of that intrinsic to the telluric correction. 

Finally, it is worth noting that from our SCUBA maps in the 
sub-mm continuum (Brand et al. 2008 ), we expect a total extinc- 
tion (i. e., over the entire width of the clump) Ay - 37 mag, or 
Ay = 18.5 mag for objects in the centre of the clump. The values 
inferred from the [Fell] line ratios towards the knots kl and k2 
are equal or less than that. It is also noteworthy that the upper 
limit to the total photospheric extinction (Ay + Ay, where Ay is 
the extinction up to the outer edge of the circumstellar envelope 
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Table 6. Visual extinction derived from [Fell] line ratios. 



Source 


Av 






1.645/1.321 


1.645/1.257 




(mag) 


(mag) 


kl 


14.0(1.1) 


14.0(0.2) 


k2 


18.5(3.2) 


14.5(3.4) 


CB230-A 


8.4(3.5) 


12.3(1.3) 



and A'y that from there to the photosphere), listed in Table IH is 
of the same order as that inferred from the sub-mm emission. 

Nevertheless, the extinction derived for the YSO spectrum, 
although within the range of total photospheric extinction {Ay + 
A'y) listed in Table |4] and obtained from the SED fits, is lower 
than 18.5 mag. This, as well as the lowest values of Ay + A'y ob- 
tained, is consistent with a picture in which the jet has removed 
part of the gas and dust from the line of sight, although we would 
expect a lower extinction from the knots, as well (see below). If 
dereddened by Ay = 8.4 mag, the intrinsic colour of the YSO is 
J-K - 2.2, which would be consistent with NIR colours of main 
sequence dwarf stars of type later than L2. Their photospheric 
temperatures are lower than the temperatures obtained from the 
SED fits (see again Table |4]l, but their spectral type would be 
consistent with the presence of water and CO absorption in the 
spectrum. However, to derive the photospheric extinction from 
the [Fell] lines we must assume that these lines originate close 
to the central protostar Instead, there are indications that other 
effects must also be taken into account. First of all, many models 
of those that describe the SED most accurately exhibit a higher 
total photospheric extinction. Secondly, the continuum appears 
redder than is usual for spectra of late-M and early-L giants, 
and spectra of dwarf stars with strong water and CO absorption 
bands (see Sect. l4.4l i. Thirdly, the extinction decreases from the 
knots to the YSO: we would however expect the YSO to be more 
embedded than the knots, since the knots are emerging from the 
cavity. And finally, once they are dereddened, the J flux density 
of knots kl and k2 is higher than that of the YSO at the shortest 
wavelengths of band J: this cannot however be true if the con- 
tinuum radiation from the knots is produced by dust scattering 
of infrared emission from the YSO. In this respect, the extinc- 
tion derived for the YSO from the 1.645/1.257 /im ratio would 
enable a more consistent picture, although we remark again that 
in the literature this ratio is found to overestimate extinction sys- 
tematically. 

In particular, the following effects may alter the intrinsic 
[Fell] line ratios: 

1. the extinction derived from the [Fell] line ratios may differ 
from the actual reddening of the continuum spectrum of the 
YSO (and the knots) because these lines might trace a differ- 
ent environment (the densest part of the jet), which is only 
seen in projection towards the YSO (i. e. matter between this 
environment and the YSO could not be accounted for in the 
extinction derivation); 

2. the observed [Fell] line ratios might be affected by absorp- 
tion lines in the background photospheric spectrum being al- 
most coincident with some of the emission lines and, then, 
lowering their measured flux; 

3 . part or all of the observed radiation from the YSO might have 
been scattered by the dust before being extincted. 

With reference to item 3, we tried to infer the effects of scat- 
tering by adopting the results of Whitney et al. ( |1997| l. Using a 
two-dimensional Monte Carlo radiative transfer code, they mod- 
elled NIR scattering in cavities, swept out by jets, of the emission 



from a protostar (surrounded by a circumstellar disk) embedded 
in an inf ailing, rotating cloud with a density profile ~ r"^. For a 
number of models, they indicated the ratios R,\ of the scattered 
flux to the intrinsic source flux {R,\ - /^s.i/^o,^) in the J,H,K, 
and L bands. We fitted a A^" function to their average values 
for Ri {Ri = 0.0013, Rh = 0.0071, and Rk = 0.018) and used 
this function to correct the spectra for scattering approximately. 
At first we assumed that the observed radiation was first scat- 
tered and then extincted; we therefore modified the intrinsic line 
ratios of [Fell] by accounting for the derived A^" "scattering" 
function. These new line ratios then replaced the intrinsic val- 
ues in the formula that defines Ay as a function of the observed 
[Fell] line ratios. This yielded negative values for Ay, showing 
either that the used set of R^ is unrepresentative of the source- 
cavity system, or that only the continuum was scattered, but the 
[Fell] emission was only extincted. 

Next, we assumed that the observed radiation was scattered 
within the cavity but not subsequently extincted. This would be 
the case if the jet had almost completely cleared the cavity and 
this was inclined towards the observer Using the same mod- 
els that were later used by Whitney et al. d 19971 1, Whitney & 
Hartmann d 19931 1 found that if the angle between the cavity axis 
and the line-of-sight is less than 66°, the image of the cavity no 
longer looks bipolar because of the extinction of its red side. In 
the case of CB230, since the cavity is not bipolar (see Fig.[rii, it 
must be inclined towards the observer, which is also suggested 
by the far smaller red lobe of the CO outflow. It is therefore 
plausible that scattered light from the YSO is little extincted in 
its path to the observer Remarkably, if we use the A'" func- 
tion instead of the reddening law and correct the spectrum only 
for scattering, we obtain a bluer continuum (see Fig. |6]l and the 
[Fell] 1.645/1.321 and 1.645/1.257 line ratios are only 7% and 
14% less than the intrinsic ones, respectively. Then, the adopted 
scattering correction would naturally yield the intrinsic [Fell] 
line ratios, suggesting that the observed YSO spectrum is also 
consistent with pure scattered emission from the innermost re- 
gions of the globule, emerging through the cavity. 

Figure |6] shows examples of the YSO spectrum after correc- 
tion for scattering only, after correction for dust extinction of 
Ay - 12.3 mag, and after correction for Ay = 8.4 mag. 

In summary, our scattering model is oversimplified. Without 
more detailed modelling of CB230-A, it is difficult to disentan- 
gle the effects of scattering and extinction. We therefore adopt 
the conservative approach of correcting the YSO spectrum only 
for a reddening of Ay = 8.4 mag, derived from the [Fell] line ra- 
tios. This choice does not affect the measured absorption-line 
equivalent widths, or the emission-line fluxes, which remain 
good estimates provided that they are produced along the jet, 
but far from the protostar-disk system. 

4.4. Absorption features towards CB230-A: a tentative 
spectral identification 

The presence of absorption features in the YSO spectrum, if pro- 
duced by the photosphere, would allow us to constrain the phys- 
ical parameters of the central protostar more reliably. Spectral 
identifications of Class I sources have been successfully carried 
out in the NIR by a number of authors (see e. g. Nisini et al. 
I2005al Doppmann et al. 120051 and references therein), but al- 
ways using medium or high spectral resolution data. These au- 
thors found ubiquitous CO overtone absorption bands, and a 
number of absorption lines, due to the "cold" photospheres of 
the central protostars. However, CB230-A is the first instance 
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Fig. 6. From top to bottom: YSO spectrum corrected only for 
scattering (see Sect. 14. 3t . corrected for Ay = 12.3 mag and cor- 
rected for Ay = 8.4 mag. All spectra have been normalised to 
the peak of the emission and shifted along the y-axis. 



of a Class I source exhibiting both water absorption bands and 
other absorption features at low spectral resolution. 

In Fig. m the spectrum towards the YSO is shown after 
dereddening by Ay - 8.4 mag, i. e. the value derived from the 
[Fell] 1.645/1.321 line ratio. The continuum emission peaks in 
the H band and there are prominent absorption bands longward 
of 1.3 fim (H2O) and 2.29 yum (first and second CO overtone). 
Other absorption lines (also labelled in Fig.|4]l were found in the 
spectrum; we identified them by comparison with spectra of late- 
M and early-L type photospheres. Both detections (S/N > 3) and 
tentative detections of absorption features are listed in Table |7] 
following the spectral catalogues of Wallace et al. ( 12000 1 MK 
standards in the 7-band), Meyer et al. ( 119981 MK standards in 
the //-band), Lancon & Rocca-Volmerange (|1992| l and Dallier 
et al. (119961 both papers discuss NIR spectra of O- to M-stars of 
luminosity classes 1 to V), Reid et al. ( 1200 11 1 and Cushing et al. 
( I2OO5. both about M- and L-dwarfs). 

Unfortunately, as can be seen in Table [T] most equivalent 
widths are only few sigma and the K-pait of the spectrum is of 
too low signal-to-noise to enable reliable detection of any ab- 
sorption lines. Nevertheless, a few lines are discernible in the 
J and H spectra. Furthermore, a few absorption features were 



Table 7. Atomic absorption lines and molecular bands observed 
towards the YSO CB230-A: detections and tentative detections 



Atom/Molecule 


A" 


EW'' 




Qim) 


(A) 


KI 


1.1778 


0.67(0.27)^ 


Mgl 


1.1828 


0.89(0.27) 


Fel, Til 


1.1974 


0.64(0.36)'^ 


KI 


1.2437 


1.76(0.29) 


KI 


1.2529'' 


0.59(0.29)'^ 


HI (Pa/3), Til 


1.2827 


1.25(0.31) 


Mnl 


1.2905 


0.70(0.31)'' 


Mgl 


1.5028 


2.63(0.62) 


OH(2-0)P(11) 


1.5575 


1.08(0.28) 


SiF 


1.5894 


1.97(0.32) 


CO(6,3) 


1.6187 


3.7(0.4) 


Fel, OH(2-0)P(15)^ 


1.6250^ 


1.4(0.4) 


CO(8,5) 


1.6617 


0.94(0.37)'^ 


OH(2-0)P(18) 


1.6877 


0.72(0.28)^^ 


Mgl 


1.7113 


1.39(0.32) 


CO(2,0) 


2.2935 


12.00(0.83) 



" All wavelengths of the atomic transitions are taken from Cushing et al. 
J2005b . except for that of Mgl at 1.18 //m, taken from the NIST atomic 
spectra database. The OH wavelengths are from Maillard et al. (1976), 
while those of CO are from Dallier et al. (fT996l| and Meyer et al. fl998] |. 

The equivalent widths were measured in the spectra after dereddening 
by using the visual extinctions derived from the [Fell] 1.645/1.321 line 
ratio (see Tab. |6l(; however, this should have little effect on the found 
values. S/N ~ 2: tentative detection. '' Blended with the [Fell] emis- 
sion line at 1.2567 yum (see Tab. [2l(. Possible contamination due to 
the OH(2-0)P(13) band at 1.5902 yum. ^ The wavelengths of Fel and 
OH(2-0)P(15) are 1.6236 and 1.6265 ^m, respectively. Here we report 
the observed wavelength. 



identified beyond any doubt: the H2O bands at 1.3-1.5, 1.75- 
2.05, and 2.3-3.2 yum and the CO(2,0), CO(3,l), and CO(6,3) 
overtone bands longward of 2.29, 2.32, and 1.62 yum, respec- 
tively. 

Given the low spectral resolution and low signal-to-noise of 
the YSO spectrum, and the uncertainty in the spectral slopes, 
we decided not to perform a spectral identification by using fits 
of template spectra. We instead used the identified features to 
check their consistency with the photospheric temperatures de- 
rived by the SED fit. These were compared with known features 
of late-type dwarfs (class V) and giant (class 111) stars: we must 
consider the latter stellar types in particular because the radius 
of the central protostar is expected to be larger, coiTesponding 
to a low surface gravity, which affects gravity-sensitive lines. 
Furthermore, the presence of a circumstellar disk causes veil- 
ing of the absorption lines, which must be taken into account. 
We also applied the diagnostics for late-type star identification 
proposed by Meyer et al. (1 19981 1 in the //-band and by Wallace 
et al. (2000) in the 7-band, and derived the values of a few spec- 
tral indices based on the water and CO bands, and Na lines, in 
late-M to L dwarfs, as devised by Reid et al. ( |2001| l, McLean et 
al. (2003 :), and AUers et al. C2007j . 

In Fig.|2l the spectrum of CB230-A (for two different extinc- 
tion and scattering corrections) is compared with low-resolution 
spectra of M3V, M8.5V, M4111, and M6111e stars. In spite of the 
uncertainties in the slopes, the H2O and CO bands ensure that 
the YSO spectrum resembles that of late M (or even early L 
dwarfs), or late M giants. Conversely, the absorption lines iden- 
tified in the J band are more typical of either early M dwarfs or 
M giants. The clearest detections are lines from KI, which are in 
fact gravity-sensitive (McGovern et al. 120041 1. The CO(6,3) band 
is also a signature of giant (luminosity class III) stars. However, 
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Fig. 7. NIR spectra of CB230-A (corrected for Ay = 8.4 
mag, upper, and for scattering only, lower) compared with low- 
resolution spectra of LHS 54 (an M3V star, from Leggett et al. 

2000) , TVLM 513-46546 (an M8.5V star, from Leggett et al. 

2001) , KV Cai- (an M4111 star, from Lancon et al. 2000) and S 
Phe (an M6IIIe star, from Lancon et al. 2000). The main lines 
and spectral features identified in the spectrum of CB230-A are 
labelled. Each spectrum has been normalised to its flux at 1.7 



the measured equivalent widths of all CO band-heads seem more 
typical of K rather than M stars (see also Origlia et al. 1993 for 
a list of typical values). These disagreements probably reflect 
both a lower gravity with respect to dwarf stars and veiling by 
the NIR emission of the circumstellar disk. Other indications of 
lower gravity come from the possible detection of the A-X 0-1 
VO band around 1.18 yum and the measured 1.14 yum Na-index 
(as defined by Allers et al. 120071 1. In summary, the YSO spectrum 
can be identified as later than M2 (with M2 a conservative limit) 
with some indications of a lower gravity than in dwarf stars. 

4.5. Cold photosphere and/or active accretion disk 

Since the mass of the protostar is expected to be < 0.5 Mq (see 
TablelU, which agrees with the mass determinations (< 0.1 Mq) 
of the circumstellar disk, clearly most of the bolometric lumi- 
nosity must be generated by the accretion. According to the 



Teff scale adopted by Meyer et al. (1 19981 1 and tabulated in their 
Table 2, Teft is ~ 3530 K for a luminosity class V (~ 3710 K for a 
luminosity class III) at the upper end of the inferred spectral type 
interval (M2). This is within the range of central temperatures 
derived from the SED fits (see Table |4]i. However, the protostar 
models of the SED-fitting tool adopt the pms tracks of Siess et al. 
(12000 ) for the central star, which, do not account for accretion. It 
would therefore be meaningful to compare our observations with 
models of accreting protostars. In this respect, Palla & Stabler 
(119911 |1992t 119931 ) follow the growth of an accreting protostar 
up to intermediate masses. Then, the evolution of the bolometric 
luminosity (including the fraction provided by accretion), Teff, 
and protostellar mass can be evaluated along their birthline. 

From our CO maps obtained at the JCMT antenna (Brand et 
al. |2008l l, we estimate an outflowing mass rate of ~ 6 x 10"^ Mq 
yr^'. Consequently, following theoretical models of star forma- 
tion (e.g. Shu et al. I2000I I. we expect an accretion rate 10-100 
times higher, i.e. around 10"^ Mq yr This is also within 
the range derived from the SED fits (see Table HJ. For such a 
value, Palla & Stabler (1993) found that a bolometric luminos- 
ity of ~ 7.7 Lq is produced by a protostar of ~ 1.5 Mq with 
TefF ~ 4467 K, i. e. a more massive central object hotter than en- 
visaged so far. We note that if we use d - 288 pc, i. e. the lower 
limit of the distance determination (Sect. [T]l, we find more pre- 
cise agreement, although still marginal, obtaining Teff - 3980 K 
and M ~ 0.6 Mq. 

Another way to improve agreement between the observa- 
tions and the model is to increase the accretion rate, which would 
produce a higher bolometric luminosity for a given mass and, as 
a consequence, T^ff could be lowered. However, if we assumed 
an accretion rate of 10""* Mq yr"' then a protostar of < 1 Mq 
would need less than 10^ yrs to be assembled. But this would 
be less than our estimated outflow dynamical age (~ 3 x 10'* yrs 
at d = 450 pc, assuming an inclination with respect to the line- 

o 

of-sight of / ~ 45 ). If the dynamical age can be considered to 
be a lower limit to the age of the driving source, this rules out 
an accretion rate of more than 10""* Mq yr This again agrees 
with the upper limit of the accretion rate derived by the SED fits. 

Wuchterl & Tscharnuter (2003) presented evolutionary 
tracks of protostars following the collapse of Bonnor-Ebert 
spheres. Their tracks showed that Tetj ~ 3000 K and bolometric 
luminosities of ~ 7.7 Lq can be obtained for central protostars 
of ~ 0.5 Mq and < 0.3 Myr old. By using these tracks, we 
therefore obtain a far closer agreement with our observations. 
Nevertheless, according to this model, the protostar would have 
already accreted ~ 90 % of the initial clump mass, which is 
inconsistent with the high current envelope mass that we have 
measured (~ 3.4 Mq). 

One other possibility to be considered is that a few spectral 
features (namely the strong absorption bands of CO and H2O) 
are not photospheric in origin, but are produced in an active ac- 
cretion disk. The NIR spectra of FU Ori stars exhibit strong ab- 
sorption bands of water and first and second overtones of CO, 
in contrast with their optical identification as F-G supergiants 
(Mould et al. 1 19781 1. This was attributed to the NIR part of FU 
Ori spectra being generated by an accretion disk (Calvet et al. 
|1991al l. Calvet et al. ( |1991bl l showed that for central stars with 
TefF < 6000 K, mass accretion rates of > 10"^ Mq yr"' are able 
to turn the CO bands into absorption. Even if the mass of the pro- 
tostar were 1 Mq, it could have been accumulated in less than 
10^ yr with an accretion rate of only > 10"^ Mq yr"' . These val- 
ues lie at the upper end of the acceptable ranges of both age and 
mass. It is therefore possible that at least some of the features in 
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the YSO NIR spectrum are produced in an active accretion disk 
rather than in a protostellar photosphere. 

Clearly, these issues will be resolved by future high- 
resolution high-sensitivity NIR spectral observations, which will 
enable a more robust spectral classification to be achieved and al- 
low us to recognise the kinematical signature of a disk in the line 
profiles, irrespective of reddening and scattering corrections. 

4.6. Spectral indicators of tlie accretion rate 

MuzeroUe et al. ( 119981 ) found a tight correlation between the 
emission luminosity of Bry and the accretion luminosity in a 
sample of classical T-Tauri stars. We can use their Eq. 2 and the 
line luminosity estimated from the flux listed in Tab. |2] (assum- 
ing d = 450 pc) to derive an accretion luminosity of ~ 0.15 L0. 
This would be consistent with "low" accretion rates (~ 10"^ Mq 
yr"') and central masses (~ 0.1 M0), but in total disagreement 
with the observed "high" bolometric luminosity, which in this 
case should be produced mostly by accretion. By raising the cen- 
tral mass to obtain higher stellar luminosities, we must however 
resort to the active disk scenario to produce the CO bands in 
absorption, which would again disagree with the low accretion 
luminosity derived from the Bry flux. Part of the discrepancy 
may be due to the likely underestimate of the extinction (or the 
neglect of scattering) already discussed. For instance, an under- 
estimate in Ay of 3 mag would increase the derived accretion 
luminosity to ~ 4.5 L0, in far closer agreement with the esti- 
mated bolometric luminosity. However, this would assume that 
the Muzerolle et al. ( 1 998 J relation is applicable in addition to 
less evolved YSOs. 

5. The outflow 

5. 1 . Pliysicai properties of tlie jet from tlie NIR spectra 

The extinction of the knots was derived in Sect. 14.31 using the 
ratios of the intense emission lines due to three [Fell] transi- 
tions at 1.257 fim, 1.321 jim, and 1.645 //m. The values ob- 
tained are listed in Table |6] and we adopted those derived from 
the 1.645/1.321 jjm ratio. All the other observed [Fell] lines (see 
Table |2] and Figs. |2] and [3) are associated with similar excita- 
tion levels (~ 11000-12000 K) and therefore cannot be used 
to derive temperature estimates (e.g. Nisini et al. 2005b). On 
the other hand, these [Fell] lines have different critical densi- 
ties ranging between lO'* and 10^ cm"^. Following the model 
developed by Nisini et al. ( |2005bl l, it is possible to derive an 
estimate of the electron density («e) by using intensity ratios in- 
volving the brightest transitions which are density-sensitive (i. 
e. 1.645/1.600, 1.645/1.533, and 1.645/1.677). For the CB230-A 
position, we obtain 1.645/1.600 > 2.2, 1.645/1.533 = 2.2, and 
1.645/1.677 = 4.3, which are consistent with «e > 10^ cm"-*. In 
other words, the density towards the YSO is higher than the crit- 
ical density of the [Fell] lines, which are therefore expected to 
be thermalised. For the knot kl, we derived 1.645/1.600 - 10.3, 
1.645/1.533 = 4.3, and 1.645/1.677 = 12.5, which correspond to 
electron densities in the 6x10^ -1x10"* cm range. Finally, the 
weaker k2 spectrum indicates 1.645/1.600 > 7.5, 1.645/1.533 = 
5.4, and 1.645/1.677 > 7.5, again suggesting densities of around 
10"* cm"''. Our measurements confirm the findings of Nisini et 
al. j2005bl l in HH-flows associated with three SFRs, i. e. that the 
NIR [Fell] emission lines trace either high density post-shock 
portions of the ionised gas or regions of a high degree of ionisa- 
tion and then a high electron density. Unfortunately, to derive the 
hydrogen ionisation fraction Xg, following the procedure intro- 



duced by Bacciotti & Eislofltel ( 119991 ), we would need measure- 
ments of oxygen and nitrogen lines in the optical and to assume 
that these lines trace the same material that emits the NIR [Fell] 
lines. To derive a rough estimate of the total hydrogen density, nu 
- njxg, we can only assume that the ionisation fraction is equal 
to the typical value measured in a number of jets (e.g. Podio et 
al. 120061 and references therein), = 0.03-0.6. If this is also 
valid in the CB230 case, we infer densities riu towards the kl-k2 
jet of between 10'* and 3 x 10'' cm"^. 

Finally, the comparison between the [Fell] 1 .64 /im and the 
H2 2.12 fim intensities indicates that these emission lines are 
not spatially correlated, in agreement with the scenario in which 
[Fell] and H2 are excited by different mechanisms, and the [Fell] 
emission originates in dissociative shocks. Furthermore, the H2 
line intensity clearly increases from the outermost knot to the 
YSO (see Table H. 

5.2. The outflow driven by the jet 

Regarding the mass-loss process, it would be instructive to com- 
pare the momentum flux of the atomic jet (Pjet) with that of 
the larger-scale CO outflow (Pqo - 10"^ Mq yr"' km s"'). 
Unfortunately, the present NIR data do not allow us to measure 
the momentum flux of the jet directly. By using our derived val- 
ues of ionisation density (i. e. in the range 6 x lO"*- 1 x lO'* cm"-') 
and jet size (FWHM = 1''2), and assuming typical values for 
the ionisation fraction of 0.03-0.6, and a jet velocity of 100- 
500 km s"', we obtained a range of values for /"jet of 10"^- 
10"^ M0 yr"' km s"'. The upper limit to this range is clearly 
unrealistic, which supports a high ionisation fraction, > 0.1, 
to limit the jet momentum flux arbitrarily to the more reasonable 
range of lO-^-lO"'* M0 yr"' km s"' (e.g. Podio et al. 200Q. 
In any case, it appears that Pjet > Pco, which indicates that the 
[Fell] jet is indeed able to drive the molecular outflow. 

6. Summary and conclusions 

We have presented our study of the star-forming site located in 
the globule CB230 using spectroscopic NIR observations. We 
obtained JH and HK spectra of the knots in the [Fell] jet, which 
originates in the low-mass YSO CB230-A, and the protostar it- 
self. This allowed us to derive the physical properties of the 
region where the jet is launched from. We have also retrieved 
Spitzer archive data in the range 3.6-70 yum that, complemented 
with our NIR spectra, enabled us to constrain the relevant phys- 
ical parameters of the protostar. 

Our main results are the following: 

1 . the spectra of the jet knots show a large number of emis- 
sion lines, including a rich set of [Fell] lines. The brightest 
[Fell] and H2 emission lines are spatially uncorrected, con- 
firming that [Fell] and H2 are excited by different mecha- 
nisms, in agreement with the models in which [Fell] traces 
dissociative J-shocks and molecular hydrogen traces slower 
C-shocks. 

2. the YSO spectrum exhibits a large number of atomic and 
molecular emission, and absorption features. The character- 
istics of this spectrum, if photospheric in origin, would cor- 
respond to a spectral type later than M2 for CB230-A. They 
also suggest a lower gravity than for class V stars, which is 
expected in a protostar where the radius is larger than that of 
a corresponding main sequence star. 

3. complementing flux measurements at different wavelengths 
with the Spitzer data, we obtained a SED that can be fitted 
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by a young low-mass (M < 0.5 Mq) protostar with a high 
( 1 0"'* - 1 0"^ M0 yr" ' ) accretion rate and a circumstellar disk 
with M < 0.05 M0. These limits agree with the constraints 
obtained from the NIR spectrum and the other available mm 
and sub-mm observations. The occurrence of a cavity of rel- 
atively small inclination with respect to the line-of-sight is 
in accordance with the NIR morphology of the object and 
explains why a spectrum of the innermost protostellar region 
has been obtained. 

4. as for CB230-B(,C), the lack of resolved FIR and sub- 
mm observations prevented strong constraints by SED fit- 
ting being achievable; a large spread in the central mass 
(0.1 - 1.5 Mq) and evolutionary stages (from Class I to Class 
II) was obtained. 

5. the absorption bands of water and CO in the YSO spectrum 
might be the signature of an active accretion disk, rather than 
an indication of a photospheric origin. NIR spectra of both 
high spectral resolution and high signal-to-noise ratio are 
needed to disentangle the photospheric emission from that 
of the circumstellar disk, further constraining the protostel- 
lar physical parameters. 

6. by using intensity ratios involving density-sensitive [Fell] 
emission lines and following the model by Nisini et al. 
(I2005bl l. we estimated the electron densities of the [Fell] 
knots along the jet, obtaining values in the range 6x10^- 
1 X 10"* cm-l 

7. the [Fell] jet appears to be able to drive the larger-scale out- 
flow, exhibiting a higher momentum flux than the CO out- 
flow. 
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Appendix A: Estimate of errors due to inaccuracies 
in the correction of the atmospheric absorption 

The eff^ect of an inaccurate correction of atmospheric absorption 
is assessed as follows. We assume that T(A, t) is the atmospheric 
transmission at wavelength A and time t. If the source spectrum 
S{A) taken at time t\ is corrected with a T estimated at time t2, 
then the "corrected" spectrum will be: 



Irrespective of the shape of T, we can then write 

T{A2,ti) TiAuh) 



T(A2,t2) T(Aut2) 



+ A 



so that Eq. |A.2| becomes: 
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The error in the line ratios should then depend only on the 
ratio of r at a given time to that at a different time. 

This ratio can be easily estimated from the ratio of (ex- 
tracted) spectra of a single object (be it the target or a telluric 
standard) taken at different times and, obviously, not yet cor- 
rected for the atmospheric absorption. We note that this ratio 
also includes variations due to large-scale sensitivity differences 
in the detector area and slit misalignments during AB cycles. 

As for our observations, we consider spectra of telluric stan- 
dard stars (HIP97033 and HIP108772), taken both 7 hours be- 
fore and just prior to the observations of the target. The different 
ratios of T depend almost linearly on A in the relevant intervals 
within the J-, H-, and /T-bands, increasing or decreasing with 
wavelength. The larger slope (in absolute value) is always found 
in the 7-band; it is much smaller for the other bands. We es- 
timated the ratios of T at 1.644, 1.32 and 1.25 yum (JH), and 
at 2.12 and 1.644 jim (HK). The quantity A derived by using 
spectra of telluric standards within a same ABBA-cycle (time 
differences of a few minutes) is almost always ~ 0.1, or even 
less. It is likely that, in this case, A is dominated by variations 
due to large-scale sensitivity differences and slit misalignments 
during the AB-cycles. When using the target spectra (time dif- 
ferences of a few to 10 minutes), we found again that A ~ 0.1. 
Finally, using telluric standard spectra separated by ~ 7 hrs, we 
obtained A ~ 0.8 - 1 for JH and A < 0.35 for HK. By roughly 
assuming that A increases linearly with time, we again found 
that A ~ 0.1 in JH (much less in HK) within ~ 1 hour, i. e. the 
maximum time interval between telluric and target observations 
in JH. Regarding atmospheric absorption, we conclude that line 
ratios for data acquired using a same grism should be accurate 
to within 10%, after taking account of both variations due to 
large-scale sensitivity differences and slit misalignments during 
AB-cycles. 

Problems arise only in observed wavelength ranges affected 
by the strong telluric absorption between the J- , H- and K- 
bands, which coiTespond to the edges of the three bands. In this 
respect, the changes in slopes occuiTing at wavelengths close to 
the edges of the three segments of spectrum (where the noise 
also increases) have to be regarded with caution. They might 
have been increased by variations of telluric absorption that have 
been unaccounted for. 
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The estimated ratio between two lines at wavelengths A2 and 
Al will then be: 



Rf.fit — 
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Appendix B: Reliability of the [Fell] 1.645/1.257 nm 
line ratio as an extinction estimator 

As stated in the text, the [Fell] line ratios 1.645/1.321 yum and 
1.645/1.257 yum can be used to estimate the extinction towards 
the [Fell] -emitting regions, since all of these lines originate in 
the same upper energy level. However, Nisini et al. ( I2005bl l 
found that the 1. 645/1. 257-ratio systematically underestimates 
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the extinction towards a number of sources. They suggested that 
this may be due to a poor determination of the Einstein coeffi- 
cients. This can be checked easily because any discrepancies in 
the adopted Einstein coefficients should produce a fixed system- 
atic difference between the values of Ay determined from the two 
line ratios for any observational dataset. If this difference was 
then found not to be constant (within errors), alternative sources 
of contamination should be considered. 

Such a check is beyond the scope of this paper We propose 
another possibility that does not involve the Einstein coefficients; 
the cause would be the same as that envisioned by Bailey et al. 
(1200 7). i. e. an error in the atmospheric correction occurring in 
low (and possibly also mid) spectral-resolution spectra, in this 
case affecting the 1.257 fim line. This line is in a wavelength 
range containing the 'Ag absorption band of O2. At high spec- 
tral resolution. Fig. 8 of Bailey et al. (2007) shows that there are 
many narrow lines that lower the transmission to almost zero 
at their central wavelengths, in the same region in which the 
1 .257 jum line is located. At low resolution, these high-opacity 
absorption lines are smoothed away and produce a broad depres- 
sion in the empirically-obtained atmospheric transmission curve. 
If the observed wavelenght range of the [Fell] line corresponded 
to one of the narrow, atmospheric absorption lines, it would be 
underestimated following telluric correction. This would lead to 
an overestimate of the 1. 645/1. 257-ratio and, in turn, to an over- 
estimate of the extinction when this was derived from the ratio, 
the situation which has indeed been found. It is difficult to assess 
the magnitude of the error without a detailed modelling of both 
the atmospheric absorption and the [Fell] line. 
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